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We study the effect of an optical lattice (OL) on the ground-state properties of one-dimensional 
ultracold bosons with three-body attraction and two-body repulsion, which are described by a cubic- 
quintic Gross-Pitaevskii equation with a periodic potential. Without the OL and with a vanishing 
two-body interaction term, soliton solutions of the Townes type are possible only at a critical value 
of the three-body interaction strength, at which an infinite degeneracy of the ground-state occurs; 
a repulsive two-body interaction makes such localized solutions unstable. We show that the OL 
opens a stability window around the critical point when the strength of the periodic potential is 
above a critical threshold. We also consider the effect of an external parabolic trap, studying how 
the stability of the solitons depends on matching between minima of the periodic potential and the 
minimum of the parabolic trap. 



I. INTRODUCTION 



In the current studies of ultracold quantum gases, a great deal of interest has been drawn to the study of Bose- 
Einstein condensates (BECs) loaded into optical lattices (OLs), i.e., spatially periodic potentials induced by the 
interference between counterpropagating laser beams P, 0) @|- Besides playing a crucial role in effectively tuning 
the interaction strength in the condensate, i.e., the ratio between the kinetic and interaction energies OLs offer 
an extremely useful tool for studies of the transition between the superfluid and Mott-insulator states 0, and for 
the investigation of effects in the matter-wave dynamics due to the interplay between nonlinearity and the quasi- 
discreteness, which is induced by a deep lattice potential Q. 

The mean-field dynamics of the BEG loaded into the OL is described by the cubic Gross-Pitaevskii equation (GPE) 
with the periodic potential P, [3, 0]- The respective Bogoliubov's excitation spectrum features a band structure, 
similar to electronic Bloch bands in solid state. If the OL potential is deep enough, the lowest-band dynamics may 
be approximated by the discrete nonlinear Schrodinger (NLS) equation Q. Using this correspondence, the BEG 
dynamics was studied in the framework of the nonlinear-lattice theory, see works 0, Q and short review 0]. 
The presence of the O L g ives raise to energetic and dynamical instabilities, which have been predicted theoretically 
[IS 111 [11, 111, [11, III, [li, [13 and studied experimentally [H, [l9| . 

An important application of the OLs is their use for the creation and stabilization of matter-wave solitons. In 
particular, the periodic potenti al g ives rise to localized gap solitons in the case of repulsive two-body interactions, 
as was predicted theoretically [20| and demonstrated experimentally [2 1| (w ith the attractive interactions, bright 
matter- wave solitons were created and observed in condensates of ''Li [22, Hs] and *^Rb [l^l atoms). More generally, 
the use of time- and space-modulated fields acting on atoms is a powerful tool for the control of soliton properties 



25 1; for instance, while the GPE without external potentials admits stable soliton solutions only in the ID geometry 
26, 1131, OL potentials can stabilize solitons in any higher dimension [5^, [H, [3l| . Unlike ID solitons, a necessary 



existence condition for their multidimensional counterparts, stabilized by means of OLs, is that the soliton's norm 
must exceed a certain threshold value. 

Another very useful tool frequently used in experiments with ultracold atomic gases is the control of the strength 
and sign of two-body interactions by means of an external magnetic field near the Feshbach resonance 1, 2]. Further, 
recent works proposed to exploit the possibility to control the strength of three-body interactions between atoms, 
independently from the control of the two-body collisions [H, [s^. One motivation for such studies is related to 
the possibility of creating new exotic strongly correlated phases in ultracold gases. Indeed, quantum phases, such 
as topological ones or spin liquids, turn out to be ground-states of the Hamiltonian including three- or multi-body- 
interaction terms, an example being fractional quantum-Hall states described by Pfaffian wave functions [s^. In a 
recent work [35^, a ID Bose gas with A/'-body attractive interactions was studied in the mean-field approximation, 
with the objective to create highly degenerate ground-states of Hamiltonians including many-body terms. For the 
three-body interactions {J\f — 3), the system is described by a quintic GPE, i.e., the respective term in the energy 
density is proportional to |'i/'J^, where ip is the single-atom mean- field wave function (in the general case, a similar 
term is proportional to |-!/)| ). 

Soliton solutions can be found for each A/", but they represent the stable ground-state, with negative energy (which 
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is defined as per Eqs. ([6]) and ([2]), see below), only for A/" = 2, being unstable excited states with positive energy at 
jV > 4. For Af = 3, soliton solutions are ID counterparts of the well-known Townes soiitons [s^, which play the role 
of the separatrix between collapsing and decaying localized states. The Townes-like soiitons with fixed norm (which is 
1, in the notation adopted below) exist only at a single critical value of the interaction strength, at which they feature 

the infinite degeneracy [35l |: all the solitonic wave functions, if'ix) = const • [a cosh (x/cr)] '^^^ , with arbitrary width a 
(see Eqs. (|12p and ([TT|) below), have zero energy but different values of the chemical potential, fi ^ . A relevant 
issue is how this infinite degeneracy is lifted by an external potential, especially by a periodic one corresponding to 
the OL ^,[3^. 

When the two-body interaction is present, the mean-field equation is the GPE with the cubic-quintic (CQ) non- 
linearity [3^ [sll. As said above, it has been shown [33] that it is possible to tune the strength of the two-body 
interactions independently from the three-body ones. In addition to that, in the framework of the effective GPE 
for the BEG loaded into a nearly ID ("cigar-shaped") trap with tight transverse confinement, an effective attractive 
quintic term appears, in the absence of any three-atom interactions, as a manifestation of the residual deviation from 
the one-dimensionality [stI. [38j. In any case, if the two-body interaction is repulsive while its three-body counterpart 
is attractive, soliton solutions to the GQ GPE can be found in an exact analytical form (in the free space), but they 
feature an unstable eigenvalue in the Bogoliubov - de Gennes spectrum of small perturbations around them [38| , while 
the instability of the Townes-like soiitons in the quintic equation is subexponential, being accounted for by a zero 
eigenvalue. 

The issue we address in this paper is the possibility to stabilize such soiitons by means of the OL potential. 
Previously the stabilization of originally unstable soiitons by means of the OL was considered, in the 2D [2I, [2^, Isoj 
and ID [36] settings alike, only for localized states of the Townes type (recently, the stabilization of 2D soiitons against 
the supercritical collapse by the OL was also demonstrated in the GQ model in 2D, with both cubic and quintic terms 
being attractive [13]). It was found that the OL with any value of its strength (i.e., with zero threshold) opens a 
stability window around the critical point corresponding to the Townes soiitons. In this work, we demonstrate that 
the OL opens a stability window for soiitons in the GQ model (with the repulsive cubic and attractive quintic terms) 
too, but only if the lattice strength exceeds a finite threshold value. 

Apart from the context of BEG, where the nonlinearity degree is related to the number of atoms simultaneously 
involved in the contact interaction, NLS equations with the power-law and CQ nonlinearities are also known as 
spatial-domain models of the light propagation in self-focusing media [39| (for a brief overview of optical models based 
on the CQ-NLS equation, including references to experimental realizations, see recent works [10, EJ). In the case of 
the cubic nonlinearity (the Kerr medium), effects of imprinted lattices on the transmission of light beams have been 
investigated both in local "^8, 42] and nonlocal [43i, [3] models. 

The paper is structured as follows. In Section II, we introduce the GQ GPE corresponding to the mean-field 
description of the ID Bose gas with two-body repulsive and three-body attractive interactions. Properties of the 
(unstable) soliton solutions to this equation are also recapitulated in Section II. In Section III, we use the variational 
approximation (VA) (see Ref. [i^ for a review) to discuss effects of the OL on the soiitons. We introduce an appro- 
priate ansatz and compute the corresponding energy. The limit of the vanishing two-body interaction is considered 
too and compared to previous results In Section IV, the stability region for the soliton solution in the presence 
of the repulsive two-body interaction and OL is determined and compared with numerical findings. The effect of an 
additional harmonic-trap potential is studied in Section V, showing that the stability region depends on the matching 
between minima of the periodic potential and the location of the minimmn of the harmonic trap. In Section VI we 
present our conclusions. 

II. THE MODEL 

The quantum many-body Hamiltonian for the ID Bose gas with M-hody contact attractive interactions is 

f +00 



where \E'(a;) is the bosonic-field operator, c > is the nonlinearity strength and 

/.o = -|^^ + yc.(x) (2) 

is the single-particle Hamiltonian, Vc^t{x) being the external pot ential. The case oiM — 2 in the homogeneous limit 
iVext = 0) corresponds to the integrable Lieb-Liniger model [iy]. For attractive interactions (c > 0), its analytical 
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solution was obtained by means of the Bethe ansatz (47| and for a large number of particles, Ntct, the energy of 
the exact ground-state solution coincides with that obtained in the mean-field approximation [48|. In the attractive 
Lieb-Liniger model, a finite ground-state energy per particle is provided by fixing product cA^tot to a constant value 
[13, El], while for A/" > 2 one has to set c{Ntot) = const 
In the Heisenberg representation, the equation of motion for field 'i>{x,t) is 



. A/'-l / . \ A/"-! 



The mean-field approximation reduces Eq. Q to the corresponding GPE with the power-law nonlincarity, 

_d^{x, t) 



ih- 



dt 



(/io-c|V(a^,t)r)v^(a;,i), 



(3) 



(4) 



where the macroscopic wave function ip(x^t) is normalized to the total number of atoms, A^tot, and the nonlincarity 
degree is related to the order of the multi-body interactions, M: 



(5) 



Thus, the usual two-body interaction (TV — 2) corresponds to a = 2, and the three-body interaction {M — 3) to 
a = A. Equation Q conserves the energy, 



E = / dxi}*{x) 



ho - 



2c 
a + 2 



mx)\ 



a+2 



(6) 



which is the classical counterpart of quantum Hamiltonian ([T]). 

In Eq. ([2]), Vcxt(a;) is the external trapping potential, which typically includes a superposition of an harmonic 
magnetic trap and periodic OL potential, Vcxt(2^) = Vnoi^) + ^Ohix), where the harmonic confining term is Vho = 



mto'^x'^ /2. We take the periodic potential as Vol 
beams which build the OL, and q = 27r/A, with A 



e sin [qx + 5) , where e is proportional to the power of the laser 



Aiascr sin (61/2); 



is the wavelength of the beams, and 9 



the angle between them (the period of the lattice is A/2). Parameter S measures a mismatch between the minimum 
of the parabolic potential (at x = 0) and the closest local minimum of the lattice potential: when 6 — {6 = Tr/2) 
a minimum (maximum) of Vol coincides with the minimum of Vno- In fact, except for Section V, we consider the 
situation without the parabolic trap (i.e., a; = 0), therefore we set (5 = in this case. 

The time-independent power-law GPE corresponding to Eq. Q is (from now on, we use normalized units, with 
h = m = 1 and iVtot = 1) 



V'(x) = Hipix), 



where fi is the chemical potential, 



Vcxt{x) 



esin^ (qx), 



(7) 



(8) 



and the norm of the wave function is 1. In the free-space cubic model (a = 2 and e = 0), Eq. ([7]) is the integrable 
NLS equation, whose multi-soliton solutions can be obtained by means of the inverse scattering method ^27||. The 
commonly known single-soliton NLS solution is 



ip{x) = A sech {kx) . 



(9) 



where fc^ = 2|/i| and yl is a real amplitude, the respective value of the chemical potential being /i — —cA? 12. For a 
general value of a > 0, the integrability is lost even in the absence of the external potential [26|; nevertheless, the 
respective single-soliton solutions can be found in an explicit form (49l . [50| . 

For the attractive three-body interactions (a — 4), Eq. ^ is the self- focusing quintic GPE, whose stationary 
version is 



'2dx2 



ip{x) = ^'ip{x). 



(10) 



For Vcy^t{x) = 0, if one fixes coefficient c in front of the interaction term, the Townes-like solitons exist for a particular 
value of the norm of the wave function (sgI . [sij . On the other hand, fixing the normalization of the wave function 
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(recall that the norm is 1 in our units) amounts, for a ^ 4, to fixing a relation among the chemical potential and the 
interaction strength so that for each c it is possible to obtain a single soliton solution (although, as mentioned 
above, these solutions provide the ground-state in the infinite system only for a < 4, i.e., for Af < 3). However, for 
a = 4 (i.e., Af = 3) chemical potential /i remains indefinite, assuming arbitrary negative values, while the soliton 
solution of the form 

^(x) = ^''^^ifl - -8a^ (11) 
\/cosh(/ex) 

satisfies the unitary normalization condition at a single (critical) value of the interaction strength (bsI . [sij , 

c = c = — . (12) 

At c = c*, all solutions pT|) share a common value of the energy, which is simply E = Q [S^H^, as follows from Eqs. 
dSl) and (Ull). 

If the two-body interaction is added to the three-body attraction, the mean-field equation is the GPE with the CQ 
nonlinearity, 



'2 



-\^+g\i^{x)\'-c\^{x)\' + V,^,{x) 



ij{x) = iii^ix). (13) 



As said above, we chiefly focus on the case of the repulsive two-body interactions, i.e., g > 0. A family of exact soliton 
solutions to Eq. with Kxt(a:) = can be obtained in the exact form [sl, [H, [sll , which, for g > 0, is 

= , , (14) 

{l+iA^)cos\Y{2^f2\^\x)-^A^ 

where ^ = g/ (4|/i|), and the maximum value of the density, at the soliton's center, is 



A simple derivation of Eq. (fT4|) is presented in Appendix A. Obviously, for 5 = solution ((T4|) reduces to Townes-like 
soliton (fTTjl . 

Imposing the above-mentioned normalization, 

-I- 00 

\tP{x)\'^dx = l, (16) 



on solution (1141). one arrives at relation 



y|tan-i (yrrw) = 1, 



(17) 



from where it follows that, for g > 0, soliton solutions with /i < satisfying normalization condition (|16p exist for 
c> c*. However, these solutions are unstable [38j (in particular, because they do not satisfy the Vakhitov-Kolokolov 
stability criterion [s^l)- In the following section we discuss how the OL can stabilize such localized solutions. 



III. VARIATIONAL APPROXIMATION 



Both for a = 2 and 4 {J\f = 2 and 3), and for the GPE with the mixed CQ nonlinearity, the presence of the 
periodic potential makes it necessary to resort to approximate methods for finding solitons. To this end, we use the 
VA (variational approximation) [i^, [s^ based on the ansatz which yields exact soliton solution (fTTjl of the quintic 
NLS equation in the absence of the external potential: 

V'a„s(a:) = , f ^ ^ . (18) 
y'cosh[x/a) 
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Here, width a is the variational parameter to be determined by the minimization of the energy, while amplitude 
A will be found from normalization condition We expect that ansatz p^ . which does not explicitly include 

the modulation of the wave function induced by the OL, may give a reasonable estimate of the soliton's energy 
for sufficiently small values of OL strength e in Eq. ([5]), cf. the known result for the 2D equation with the cubic 
nonlinearity (a — 2) and OL potential [28,,^,]. In the case of the 3D GPE which includes the cubic term and harmonic 
trap, this approach leads to an estimate for the critical value of the number of atoms above which the condensate 
collapses, that was found to be in a reasonable agreement with results produced by the numerical solution of the GPE 
[sgI . [syj . In ID, the VA based on the Gaussian ansatz also provides for quite an accurate approximation to exact 
soliton solution ([9] [58j. Similar analyses carried out in the ID model including the cubic term and OL 0, [s^ 
have demonstrated that (unlike the 2D and 3D cases) the ID soliton trapped in the OL potential does not have an 
existence threshold in terms of its norm (number of atoms). 

The energy to be minimized in the framework of the VA is obtained by inserting ansatz (jlSp in the GPE energy 
functional given by Eq. ([6]). The kinetic and quintic-interaction energy terms in the functional both scale as cr"^; 
then, the energy per particle computed from expression ([6]) is 

i?=4 + ^ + ^[l- sech(7rga)] , (19) 



where c* is defined in Eq. (|T2|) . 

For e = (without the OL), the scenario discussed in the previous section for the uniform CQ GPE with the 
attractive three-body and repulsive two-body interactions is recovered, as energy (jl9p reduces in that case to 

+ (21) 

For 5 = 0, the energy is positive when c < c* (i.e., /3 > 0, see Eq. (f20l) ') and vanishes at ct ^ oo; for c = c* (i.e., 
/3 = 0) one obtains i? = 0, in agreement with the above-mentioned exact result showing the infinite degeneracy of 
soliton family (|lip . while for c > c* the energy is negative and diverges (to — oo) at cr ^ 0, signaling, in terms of the 
VA, the onset of the collapse. With g > 0, expression (pij) does not give rise to any minimum of the energy, which 
agrees with the known fact of the instability of all the solitons in this case [1^ . 

A detailed study of minima of variational energy (|19p is presented in Appendix B. In the following subsection, we 
consider the case of the self- focusing quintic GPE in the presence of the OL (e > 0, g — 0), while the discussion of 
the general case (e > 0, <? > 0) is given in Section IV. 



A. Self-focusing quintic GPE with the optical-lattice potential 

Here we address the stability of localized variational mode (fT8|) . for different values the OL parameters, strength e 
and wavenumber q, keeping g = 0. The results of the analysis of minima of the variational energy (jl9p . presented in 
Appendix B, can be summarized as follows (see also Fig. [T]): for c > c*, the infinitely deep minimum of the energy is 
obtained at cr ^ 0, which corresponds to the collapse, as shown in Fig. [TJa). For c < c*, the collapse may be avoided, 
and three possibilities arise: there exists another special value, c' < c*, such that for every c between c' and c* the 
energy has a minimum at a — a i and a maximum at cr = (T2, while for c < c' the energy does not have a minimum 
at any finite value of cr, see Fig. [IJd). Actually, two different situations should be distinguished for c' < c < c*: there 
exists a specific value (refer to Appendix B), 

c** = c* - -^Tc, where Tc « 2.13, (22) 
zq'' 

(with c** > c') such that, for c** < c < c*, the energy has a global minimum at cr = cti (which, thus, represents the 
ground-state of the boson gas in this situation), while, for c' < c < c**, the energy minimum at ct = tri is a local one. 
In other words, taking into regard the fact that, as shown by Eq. p^ . the energy-per-particle approaches value e/2 at 
large tr, we conclude that, for c** < c < c* (c' < c < c** ), the energy satisfies inequality E{ai) < e/2 (-E(cri) > e/2), 
as showed in Figs. [Ub,c). 

From the above analysis, we infer that for c < c** the ground-state is a delocalized one (although the metastable 
state, corresponding to the above-mentioned local energy minimum, exists for c' < c < c**), for c** < c < c* the 
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FIG. 1; Variational energy E obtained in the framework of the quintic GPE versus a (in units of e/2) for c > c* (a); c** < c < c* 
(b); c' < c < c** (c); c < c' (d). In (a) the sohd (dotted) hne is the energy for c > c* (c — c*); in (b)-(c), points of the energy 
minimum and maximum, ai and (T2, are indicated. 




FIG. 2: The numericaUy found ground-state of the quintic GPE with the periodic OL potential, for several values of nonlinearity 
coefficient c. Solid lines, starting from the narrowest configuration, refer to c = 3.7, 3.5, 3, 2.5, 2, 1.9 (recall that c* — 3tt^ /8 ~ 
3.701), and the dashed line refers to c = 1.8. Parameters are e = 6, q = 3 and L = 10. Inset: squared width of the ground- 
state as a function of c (the dot-dashed line is a guide to the eye). Critical value c** obtained from the numerical analysis 
is c** — 1.87(3), which should be compared with the corresponding value (|22l) predicted by the variational approximation, 
c** ~ 1.57. 

ground-state is represented by a finite-size soliton configuration (in agreement with Ref. [s^) and for c > c* it is 
collapsing. Equation (22) shows that the width of the stability region depends on ratio e/q^: keeping fixed all other 
parameters, the decrease of the lattice spacing (i.e., the increase of q) leads to a reduction of the stability region. 
Equation (1^ also shows that for e/g^ = 2c*/3Tc « 1.16 the VA formally predicts c** = 0: however, for c = 0, 
the ground-state is delocalized and the variational ansatz (fT5)) cannot be used, as it does not take into account the 
modulation induced by the deep OL potential. 

In Fig. [2l we plot the numerically found ground-state of the quintic GPE in a ID box {—L < x < L). It is seen 
that, with the increase of c* — c, the configuration becomes broader, until a critical value is reached, as discussed in 
[36| . In the inset of Fig. Owe plot the squared width cr^ = dx x^ \ ip{x) p of the numerically found ground-state 

versus c, which makes the delocalization transition evident: for c < c** the width cr is oc i, while around c ~ c** 
the width suddenly decreases. Variational estimate (|22|) for the critical value c**, as predicted by the VA (see Eq. 
(1^^ ). is displayed in Fig. O together with numerical results. One observes a reasonable agreement between them, 
especially for small e, which is due both to the use of the more adequate ansatz (|18p . rather than a Gaussian, and 
also because is found as the value at which the global (rather than local) minimum disappears. 
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FIG. 3: The dotted line: the variational estimate for c** as a function of e/g^, according Eq. (|22[) (for the quintic GPE), 
the dashed line corresponding to c* = Stt^/S. Discrete symbols represent results obtained from the numerical solution of the 
quintic GPE. They designate the transition form the localized ground-state to the extended one (parameters are the same as 
in Fig. ((2}. According to the variational approximation, the ground-state is delocalized (ct oo) below the dotted line, and it 
collapses (cr 0) for c above the dashed line. 



IV. THE STABILITY REGION FOR THE CONDENSATE WITH COMPETING TWO- AND 

THREE-BODY INTERACTIONS 

The most interesting situation occurs when the two-body repulsive interaction (17 > 0) competes with the attractive 
three-body coUisions (c > 0). As said above, all solitons in the free space (e = 0) are strongly unstable in this 
situation ^8!| , and the possibility of their stabilization by the OL was not studied before. The analysis of variational 
energy (|19p. presented in Appendix B, yields the following results for this case. For c > c*, the energy does not have 
a minimum at finite cr, hence the OL cannot stabilize the solitons in this case. If c = c*, the energy has a global 
minimum at a finite value of cr, when 

G = ^ < ~ 0.663. (23) 
7re 

For c < c*, the energy features a global minimum at finite cr for c**(G) < c < c*, where the modified critical value is 



c**(G)^c*-^r,(G), (24) 



cf. definition ((22)l for G = 0. The value Tc depends upon G, vanishing for G larger than the critical value Gc- This 
means that, to balance the destabilizing effect of the repulsive two-body interactions, the strength of the periodic 
potential, e, must exceed its own critical value, 

ecrit = -pr- (25) 

Otherwise, Eq. ([M)) yields c** = c*, i.e., the OL cannot stabilize the solitons. 

In Fig. [3] we plot the numerically found ground-state of CQ GPE for several values of e. It is seen that, at 
small e, the wave function ip remains delocalized, until a critical value is reached. In the inset of Fig. [4] the squared 
width of the numerically generated ground-state is plotted versus e. In Fig. [5l we compare critical value ecrit, as given 
by Eq. (|25|). with numerical results: for small g, the predicted linear dependence of Ccrit on g is well corroborated by 
the numerical results, the relative error in the slope being ~ 20%. In principle, the comparison between variational 
estimate (1251) and numerical results might be further improved by choosing a variational wave function which, in the 
limit of e = (uniform space) would reproduce exact CQ soliton (|14p . However, the calculations with such an ansatz 
are extremely cumbersome. 
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FIG. 4: The numerically found ground-state of the cubic-quintic GPE for several values of e. Solid lines, starting from the 
narrowest wave function, refer to e = 4.0, 3.0, 2.8, 2.6, 2.4, and the dashed line refers to e = 2.3. Parameters are c = 3.65, 
g = 1, q = 3, L = 5. Inset: the squared width of the ground-state versus e (the dot-dashed line is a guide to the eye). Critical 
value Ecrit obtained from the numerical data is ecrit ~ 2.35(5), which should be compared to the variational prediction given by 
Eq. which is tcrit ^ 2.88. 




FIG. 5: Solid line: ecrit versus g, as given by Eq. (|25[) . Symbols refer to results obtained from the numerical solution for the 
ground-state of the cubic-quintic GPE. They represent the delocalization transition. The parameters are the same as in Fig. 
©• 

V. THE EFFECT OF THE HARMONIC TRAP 

In this section we aim to use the variational approximation based on ansatz p8)) for examining the combined effect 
of the parabolic trapping potential acting along with an OL, i.e., we take Eq. (jlOp with external potential 

Vo^t{x)^oj^x^/2 + esin'^{qx + S), (26) 

cf. Eq. ([8]), and disregard binary collisions {g = 0). Value S = {S — tt) corresponds to the matching (largest 
mismatch) between the minimum of the harmonic potential and a local minimum of the lattice potential. The 
respective variational energy is obtained from ^ with potential (j26p : 

S = 4- + ^ ^ °" + - [1 - cos (2(5)sech(7r(jo-)] . (27) 
(J" 8 2 

With cos(25) > 0, the soliton is stable for c < c*, and it collapses otherwise. With cos(25) < 0, a richer behavior is 
predicted by the VA. The system does stabilize for c < c*, while, for c > c*, the presence of the mismatched harmonic 
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FIG. 6; The critical line separating in the {q,c — c*) plane of the model (including the parabolic trap) the metastable region 
from the unstable one. The parameters are e — 1, cos {2S) — —0.5, and u — 1. 



trap gives rise to a metastability region. Since E — > — oo as a and E — > +oo as cr — + oo, one can encounter 
two possibilities: either dE/da is positive for all a (and there are no energy minima), or equation dE/da — has 
two roots, corresponding to a local minimum and a maximum. The equation for the value of a at which energy (|27p 
reaches the local minimum is 

i/^i = '-^^^m. (28) 

where = Trqcr, and 

/ sinh 9 



q3 



Vcosh^6' 



-VO], (29) 



- 2eg2 I cos{26) \ ' 

One can see that, for c = c* (i.e., /3 = 0), Eq. (|28|) does not have a nonvanishing solution if q is smaller than a critical 
value, 

(„) ^ ^ (31) 

V2e|cos(25)| 

while it has a nonvanishing solution for g > q^"'\ 

Actually, for c > c* (i.e., P < 0), Eq. ([28|l with q > q'^"'' has two nonvanishing roots, one of which is a local 
minimum, while such roots do not exist for q < q^"''> . For q > q''"'\ the right-hand side of Eq. (f28|) has a maximum 
value, which fixes the maximum value of /3, i.e., the maximum value of c, which we refer to as c***. Then, for c > c***, 
the variational energy does not have a local minimum. For c* < c < c*** there appears a finite metastability region, 
in terms of wavenumber q, as illustrated by Fig. [6l In other words, for fixed c, metastable states appear at large 
values of e. 



VI. CONCLUSIONS 



In this work we have studied the effect of the OL (optical lattice) on the ID Bose gas with attractive three-body 
and repulsive two-body interactions, described by the GPE (Gross-Pitaevskii equation) with the CQ (cubic-quintic) 
nonlinearity. Actually, the effective quintic attractive term in the GPE may be induced by the residual deviation of 
the condensate, tightly tra ppe d in a cigar-shaped confining potential, from the one-dimensionality (when the three 
bod^^sses are negligible) |37l . [sH] or by three-body interaction terms between atoms according to recent proposals 
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FIG. 7: Plot of potential V{ipo) for ^ < 0. 



In the absence of an external potential, soliton solutions to this equation with the CQ nonlinearity are known in 
the exact form, but they all are strongly unstable. We have demonstrated that the OL opens a stability window 
for the solitons, provided that the OL strength, e, exceeds a finite minimum value. The size of the stability window 
depends on e/q^, where q is the OL's wavenumber. We have also considered effects of the additional harmonic trap, 
finding that, if the quintic nonlinearity is strong enough (c > c*), a metastability region may arise, depending on the 
mismatch between minima of the periodic potential and harmonic trap. 



APPENDIX A: LOCALIZED SOLUTIONS OF THE CUBIC-QUINTIC GROSS-PITAEVSKII EQUATION 

Assuming that ip{x) is real, we look for localized solutions to the CQ NLS equation, 

-^0 + #'-c^' = mV' (A1) 

with c > and g > 0. Interpreting a; as a formal time variable and tp(x) as the coordinate of a particle, Eq. (jAip 
formally corresponds to the Newton's equation of motion of this particle, 

where the effective mass is M = 1/2, and the potential is 

nV') = ^^'- V + (A3) 
2 4 5 

with an arbitrary additive constant chosen so as to have T^(0) = 0. Potential (|A3[) for /i < 0, which corresponds to 
normalizable solutions, is plotted in Fig. [71 Condition V{zkA) = yields expression for the soliton's amplitude. 
Further, we make use of the conservation of the corresponding Hamiltonian, 

The boundary conditions for localized solutions, ip{x oo) 0, dilj/dx{x oo) 0, select H — Q in Eq. (jA4p . 
Taking into regard the fact that < for < i/) < yl, and looking for solutions with d-^/da; < at x > 0, one 

obtains from here the soliton solution in an implicit form, 

(A5) 



11 



It further follows from Eq. (|A5p that 

with £ = e-2^2|p|2;^ -gq^ (|A6p . we use notatfon = 3|^|/c and 6^ = 3g/2c. Thus, from Eq. (|A6p one obtains 

[2a2 + 52^2 (1 _ £)]2 _^ 4a2^4£2 ■ 

One can easily check that this expression yields 'ilP'{x) = / cosh {2yj2\^\x) for 5 = 0, and that ?/'(0) = ^, as it must 
be. Finally, using relation a? + b'^A^ = A'^, one obtains Eq. from Eq. (|A7p . after a straightforward algebra. 



APPENDIX B: THE VARIATIONAL ENERGY 



In this Appendix we aim to study minima of variational energy (|19p . When g — 0, one sees that, for c > c*, 
the energy per particle tends to —00 at ^ 0, and to e/2 at cr 00. Then, with regard to dE/da > 0, no local 
(metastable) minima exist, and variational wave function (|18p is not the ground-state for any finite width. For c = c*, 
one obtains the global minimum at cr = 0, which implies the collapse. For c < c*, the situation is different: E 00 
as CT — > (because /3 > 0), and £' — e/2 — > +0 for ct — > 00. Then, it is necessary to find the value of (3 at which 
derivative dE/da has two real zeros. Introducing the parameter 

T=——^, (Bl) 

with (3 defined as per Eq. ((20)l . one can write condition dE/da = as 

„ sinh0 , 

T^e^ ^, B2) 

cosh^e ^ ' 

where 6 = trqa, as defined above. Equation (jB2p can be satisfied if T is smaller than a maximum value, T' « 2.67, 
and it then has two roots, 9i and 62, which correspond, respectively to the minimum at ct = cti, and maximum at 
CT = CT2 (see Fig. [T]). For T > T', Eq. (jBip has no roots, hence the variational energy has no minima at finite values 
of the soliton's width, a. A plot of 9i as a function of T is presented in Fig. [HI where the maximum value of 9i is 
^max ^ 3.0415. The energy minimum at 61 is a global one if E{9i) < e/2; using Eq. (fT9| . this condition reads 

- ^^'^^^ <0. (B3) 



cosh6li(r) 



As one can see from Fig. [51 condition (|B3p is satisfied for T < Tc, where Tc « 2.1289; then, a global minimum exists 
only for < T < Tc, while for Tc < T < T' the minimum is local, corresponding to a metastable state. Using the 
value of Tc and definition (jBip , one arrives at Eq. (|22p . 

For g > (recall it corresponds to the two-body repulsion), variational energy for c > c* does not have a 
minimum at finite values of a. However, for c — c* a finite minimum is possible. Indeed, with definition of G as per 
Eq. ([23)) . condition dE/da = can be written as 

G = 9'^. (B4) 

cosh^ e 

For G < G' 1.0341, Eq. (|B4p has two roots. By imposing the condition that the value of the energy a.i a = ai 
be smaller than e/2, one gets G < Gc — 0.6627. Then, similar to the situation considered above, a global minimum 
exists only for < G < Gc, while for Gc < G < G' the minimum is local. 
For c < c*, condition dE/da = reads 



T = 9'^-G9. (B5) 
cosh^ 9 ^ ' 



One can see that condition ([B5| is satisfied for T < T'{G), with T'(G') = 0. Then, for G > G', i.e., for e small 
enough, the variational energy does not have a minimum. Imposing the condition that the minimum is global leads 
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FIG. 8: The solid line represents Oi as a function of parameter T (defined in Eq. (IB2|l l for g = . The dashed line is the plot 
of function T - 261? (T)/ cosh 6*1 (T) versus T. The maximum value of 6*1 at T = T' is indicated. 




G 



FIG. 9: The solid (dashed) line is the plot of Tc (T') as a function of parameter G. 

to T < Tc, with Tc{Gc) = 0. Then, for G > Gc, i.e. for e smaller than a critical value, the variational energy cannot 
have a global minimum at a finite value of a, i.e., localized states cannot realize a global minimum. Functions T'{G) 
and Tc(G) are plotted in Fig. El in Fig. [lOl we plot maximum value 0"^^^ of 9i for T — T'{G), as a function of G. 
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